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Within a mean-field approximation, the ground state and finite temperature phase diagrams of the two- 
dimensional Kondo lattice model have been carefully studied as functions of the Kondo coupling J and the 
conduction electron concentration ric. In addition to the conventional hybridization between local moments and 
itinerant electrons, a staggered hybridization is proposed to characterize the interplay between the antiferromag¬ 
netism and the Kondo screening effect. As a result, a heavy fermion antiferromagnetic phase is obtained and 
separated from the pure antiferromagnetic ordered phase by a first-order Lifshitz phase transition, while a con¬ 
tinuous phase transition exists between the heavy fermion antiferromagnetic phase and the Kondo paramagnetic 
phase. We have developed a efficient theory to calculate these phase boundaries. As tic decreases from the 
half-filling, the region of the heavy fermion antiferromagnetic phase shrinks and finally disappears at a critical 
point n* = 0.8228, leaving a first-order critical line between the pure antiferromagnetic phase and the Kondo 
paramagnetic phase for ric < n*. At half-filling limit, a finite temperature phase diagram is also determined 
on the Kondo coupling and temperature (J-T) plane. Notably, as the temperature is increased, the region of the 
heavy fermion antiferromagnetic phase is reduced continuously, and finally converges to a single point, together 
with the pure antiferromagnetic phase and the Kondo paramagnetic phase. The phase diagrams with such triple 
point may account for the observed phase transitions in related heavy fermion materials. 

PACS numbers: 75.30.Mb. 71.10.Hf, 71.30.+h, 75.50.Ee 


I. INTRODUCTION 

Since the discovery of heavy-fermion materials, the un¬ 
derlying mechanism controlling these rare earth or actinide- 
based compounds has continuously been the focuses of ex¬ 
ploration [1-3]. In these materials, the strong coupling limit 
of the on-site Kondo spin exchange causes the Kondo screen¬ 
ing (KS) of local moments by the conduction electrons, yield¬ 
ing a Kondo paramagnetic (KP) phase. On the other hand, in 
the weak coupling limit, the Kondo coupling generates an in¬ 
direct Ruderman-Kittel-Kasuya-Yosida interaction among the 
local moments, resulting in either antiferromagnetism (AFM) 
around the half-filling of the conduction electrons, or ferro¬ 
magnetism (FM) far away from half-filling [4, 5] . 

However, in the intermediate Kondo coupling region, the 
competition between KS and magnetic correlation may pro¬ 
duce a coexisting (CE) phase with AFM and KS near half¬ 
filling [6-10]. The CE phase or so-called heavy fermion 
antiferromagnetic phase (HEAEM) has been observed in 
CeCoGe 3 _ 2 ;Sia; [1 l],and Ce 3 Pd 2 oSi 6 [9], etc. In Ce 3 Pd 2 oSi 6 , 
the HEAEM is observed within the magnetically ordered 
phase, indicating the separation of two transitions, i.e., 
the breakdown of Kondo screening effect and concomitant 
Eermi surface reconstruction (ESR) which happens between 
HEAEM and pure AEM, and the magnetic transition which oc¬ 
curs between HEAEM and KP phase [9]. However, studies of 
Hall coefficient and pressure effect in YbRh 2 Si 2 have shown 
that the Kondo breakdown occurs precisely at the magnetic 
transition, while under Co and Ir doping, two transitions sep¬ 
arate [7, 8, 10]. The Kondo breakdown also occurs away from 


the magnetic transition in Celns and CeRhi_a;Co 3 ;In 5 [12, 
13]. 

To understand the novel behavior of the phase transitions 
in these materials, the corresponding parameter region and 
the feature of the AEM phase, CE phase, KP phase have 
been investigated within the framework of the Kondo lattice 
model (KLM) or Kondo Heisenberg lattice model, and in¬ 
tensively studied at zero temperature by mean-field approx¬ 
imation, variational Monte Carlo calculations, Gutzwiller ap¬ 
proximation, etc [6, 14-18]. At half-filling limit, the recon¬ 
struction of the energy bands leads to an insulating state, and 
the ground state phases were computed with varying Kondo 
coupling [6, 14, 15]. Away from half-filling, the phase transi¬ 
tions are discussed by mean-field approximation on the Kondo 
Heisenberg lattice model, and the shift from onset to offset 
between the Kondo breakdown and magnetic transition is pro¬ 
posed to be driven by the change of the Heisenberg exchange 
and the ratio of short and long-range hopping strength [10]. 

Actually, the ground state phases and their features are con¬ 
trolled by the Kondo coupling J, the electron occupy num¬ 
ber Tic and the electron hopping strength. However, the phase 
evolutions with these parameters have not been fully explored 
yet, particularly how these phases evolute with ric and long¬ 
distance electron hoping remains an open issue [5, 18]. In 
order to deal with this issue, we adopt the slave-fermion 
mean-field technique on the KLM, and developed a more effi¬ 
cient theory to calculate the phase boundaries between various 
phases. We show that the relative positions of Kondo break¬ 
down and magnetic transition can be shifted by both ric and 
t' on the ric-J plane. As ric decreases, two transitions get 
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closer and then coincide. This ric-generating offset-to-onset 
structure with a triple point in the phase diagram is related to 
the experimental observations of Ce 3 Pd 2 oSi 6 and YbRh 2 Si 2 
under doping. 

On the other hand, most of earlier works focused on the 
ground state, while detailed theoretical studies at finite tem¬ 
peratures are still lacking. In the case when the occupation 
number of the conduction electron ric is far away from half¬ 
filling, the coexisting phase of FM and KS has investigated 
at finite temperatures and the phase diagram has been de¬ 
rived [19, 20]. Remarkably, the boundaries separating the 
various phases joint to a single point on the J-T plane, so 
an interesting issue arises as to whether the half-filled Kondo 
lattice system possesses similar feature in the finite tempera¬ 
ture phase diagram. Therefore, we devote to study the ground 
state and finite-temperature phase diagram of the Kondo lat¬ 
tice model at and away from half-filling, attempting to give a 
detailed description of the evolution of the various phases with 
the Kondo coupling, the conduction electron occupy number, 
the electron hopping strength and temperature. To this end, a 
modified mean-field decoupling technique for the Kondo in¬ 
teraction is employed, then the phase diagrams are obtained 
as functions of ric, J, t', and T . Our method turns out to give 
a compact description of the phase diagrams at both zero and 
finite temperatures. 


II. MODEL AND MEAN-FIELD TREATMENT 

We consider the spin-1/2 Kondo lattice model on a two- 
dimensional square lattice with N sites, 

H = ^(Ck - fijCkcrCkcr + S, • Sic, (1) 

k.cr i 

where ek is the dispersion of conduction electrons, which in¬ 
teract with local moments through antiferromagnetic Kondo 
exchange J > 0, and p denotes the chemical potential. 
Sic = 5 I]a /3 cL'^«/ 3 Ci /3 with (T being the Pauli matrix, rep¬ 
resents the spin density for conduction electrons, while the 
local moments can be written in the slave-fermion represen¬ 
tation as Si = i Yjap fL^apfip, which is subject to the 
restriction fiaficr = 1 imposing by a Lagrangian term 
Si ~ !)■ The Kondo interaction can be de¬ 

composed into 

S,. s,c = - g 

i i 

+ I - fy^l) 

i 

+ I J2^4tfii + 4lfit)ifyit + /ifC^f) 

i 

+ I 

i 

where the first term represents the Kondo singlet screening 
effect and the other three terms describe the triplet parings 


between conduction electrons and slave fermion holes. This 
expression captures SU(2) invariance of the Kondo coupling. 

In order to describe the antiferromagnetism in the Kondo 
lattice model, two AFM order parameters 

<7 <7 

( 2 ) 

are introduced to decouple the longitudinal Kondo spin ex¬ 
change coupling [6], where Q = (tTjTt) is the AFM vec¬ 
tor. Then the total staggered magnetization is expressed by 
M = nif — rric- To characterize the KS in the presence of 
AFM long-range ordering, two different hybridization param¬ 
eters on each magnetic sublattice A and B have to be intro¬ 
duced [6, 18] 

The conventional hybridization parameter is expressed as 


K = ^iVi+V2) 

~ 4AlfiAi) — 2^4B^fiB^ + 4BifiBi), 

while the staggered hybridization parameter is defined by 


V, = 


2(^1 - V2) 


-1 




which requires the breaking of particle-hole symmetry of the 
conduction electrons, i.e., ric ^ 1, or /' ^ 0. It is seen that 
the singlet channel hybridizes the c- and /-fermions with the 
same wave vector, while the longitudinal exchange brings a 
momentum transfer Q within both c- and /-fermions, result¬ 
ing in the staggered triplet channel. The local Lagrangian con¬ 
straint is replaced by a uniform one: \i = A. 

Though such a mean-field treatment, the model Hamilto¬ 
nian is written in the momentum space by the matrix form 

7/ = iVeo + X'^L^ka^ka, (3) 

k,fT 


where the superscript represents the summation of k restricted 
in the magnetic Brillouin zone (MBZ) with boundaries | ± 

ky\ = TT, a four-component Nambu operator has been used 
T'ko- = (ck<T Ck+Q<T ha h+QaV, and the constant term is 
given by eo = -^ (31/,^ ~ + J'mc'm / — A. The Hamiltonian 
matrix is given by 

/ ek-/i ^Jmfa -\JVs jJVta \ 

_ ^Jrrifa Ck+q -/i jJVta -\JVs 

-iJVs \JVta A -\Jm,a ’ 

\ \JVta -jJVs -^JuicCr A / 

where Ck = —2t{coskx + cosky) + At' cos k^ cosky is the 
tight-binding dispersion of conduction electrons with nearest- 
neighbor (NN) and next-nearest-neighbor (NNN) hoping 
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strength t and t', respectively. In general filling case of con¬ 
duction electrons, the quasiparticle excitation spectrums can 
not be derived analytically, we have to perform numerical cal¬ 
culations. However, at half-filling, the particle-hole symmetry 
can help to simplify the related calculations. 

III. ZERO-TEMPERATURE PHASE DIAGRAM AT 
HALE-EILLING 

We first discuss the half-filling limit with only NN hop¬ 
ing t. In this case, the particle-hole symmetry guarantees 
X = n = 0. Moreover, the staggered hybridization disap¬ 
pears as Vt = 0. Further discussions including the influence 
of Uc and NNN hoping t' will be given in the last section. 
The NN hoping between conducting electrons leads to the dis¬ 
persion ek = —2t{coskx + cosky), satisfying the relation 
Ek-fQ = —Ek- For convenience, we define V = —§14, then 
the analytic formulas for the four dispersions are obtained by 
diagonalizing 'Hka- 



FIG. 1: (Color online) (a) (b) Zero-temperature energies of AFM 
phase KP phase and CE phase E^^ vs Kondo coupling 
J. (c) Free energies for the three phases as functions of J at temper¬ 
ature T = 0.2f. The CE phase exists in a narrow range between Jd 
and Jc 2 . All energies are in unit of NN hopping strength t. 


where 





phase corresponds to the phase with lowest energy. When 
rric = mf = 0, the energy density of KP phase is obtained 
by 



E2k = y-F ^ J4(mcm/ -F V^)^, 


with the relation 

F'k ~ \/F'lk + £'2k = F'k- 

At zero temperature, two lower branches of spectrums 
—E^ lying below the Fermi level are full occupied, giving 
rise to an insulating heavy fermion state. Two higher branches 
E^ above the Fermi level give no contribution to the ground 
state energy, therefore the ground state energy density is eval¬ 
uated as 


(5) 


with the summation of k runs over the entire Brillouin zone. 
The mean-field order parameters are determined by minimiz¬ 
ing E^^, yielding the self-consistent equations 


J 

I/V 

8N 

3J 


T— = 


i?k 3(m/ - me) 


{mcmf +V^)Y^ 


EkE2k 


E 


£^kF^2k 


= 1 -F mf /me- 


2m f — 3mc 
3{mf — me) ’ 


( 6 ) 


In the whole coupling range, the pure AFM phase and KP 
phase should also be examined, then the stable ground state 


E. 


•KP 



1 

N 




(7) 


where V is determined by 4-V. —= 4/3. For 

J' N ^k ^e2^j2y2 / 

the pure AFM phase with U = 0, the conduction electrons 
and the local moments are totally decoupled, resulting in the 

dispersions E^ = ± -F for the conduction elec¬ 

trons and E^ = ±^Jme for local moments, respectively. By 
performing similar self-consistent treatments, the energy for 
AFM phase is found to be 


^E = -^EV^k + ^vi6, (8) 

k 

with order parameters mf = 1I2, and me = 

J_ V 1 

The comparison of the ground state energies for AFM, CE 
and KP phases is demonstrated in Fig. l(a)-(b). As expected, 
the competition between ordering of local moments and for¬ 
mation of Kondo singlets leads to a coexisting solution with 
lowest energy, indicating the stability of the CE phase in the 
intermediate Kondo coupling range < J < Jc 2 [6], while 
the pure AEM phase and the KP phase exist in the region 
J < Jci and J > Je 2 , respectively. The derived staggered 
magnetization M = mf — me and KS strength V are given 
in Pig. 2(a) as a function of the Kondo coupling J. In the 
CE phase, the fluctuations of local spins in the Kondo channel 
suppress the AEM order, while the staggered magnetic order 
brings a rapid decrease of the KS strength. Both M and V 
vary continuously on the phase boundaries. and 4^2 cor¬ 
respond to the lower boundary of the KS state with order pa¬ 
rameters me 0,mf 0,V —>0, and the upper boundary of 
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FIG. 2: (Color online) Staggered magnetization M and Kondo hy¬ 
bridization V as functions of Kondo coupling J at (a) T =0 and (b) 
T=0.2t for the three phases. 


the AFM state with rric —>■ 0, my —>■ 0,1^ 0, respectively. 

The limit 1^ —s> 0 can be replaced by setting 1/ = 0 in the self- 
consistent equations Eq. (6), because the denominators of the 
integral functions are always nonzero, leading to the numeri¬ 
cal results Jci = 3.1498f, my = 0.5, rric = 0.2735. To calcu¬ 
late Jc 2 , we can simply set rric = 0, my = 0 in the integrals 
in Eq. (6) because E’k = \/ek + >^ 02 ^^ E 2 k = 
are both gapped. To keep my /me as a constant, we obtain the 
numerical solutions Jc 2 = 3.4161f, V = 0.5853. 

As noticed, the energy of CE phase is tangent to 

E^^ on the edge Jd, and to E^^ at Jc 2 , respectively, im¬ 
plying that the Kondo lattice system undergoes second-order 
phase transitions on both phase boundaries. The coexistence 
of AEM and KS in the Kondo lattice systems has been re¬ 
ported previously at zero temperature [6, 10, 14]. Though the 
proposed phase boundaries Jd and Jc 2 are slightly different 
from our results, the main physical pictures remain the same. 
This inconsistence comes from distinct mean-field decoupling 
procedures. In these earlier studies, the Kondo exchange is de¬ 
composed directly into longitudinal term and transversal part 
(which is the singlet channel of hybridization), respectively; 
while in this work, the singlet and triplet hybridization be¬ 
tween itinerant electrons and local moments are considered at 
the same level in the beginning. Therefore, our method can be 
generalized straightforwardly to deal with the case away from 
half-filling and with hoping beyond NN, as will be discussed 
in the following. 


IV. FINITE TEMPERATURE PHASE DIAGRAM AT 
HALE-FILLING 


the free energy density can be calculated via the partition 
function, leading to the result 

2 2r 

= -JE^+Jmemy-—^ln[2cosh(j;^/2r)]. (9) 

k.± 


It is easy to verify the equivalence of the ground state energy 
and above free energy at zero-temperature limit. The mean- 
field parameters are determined by minimizing F^^, then the 
self-consistent equations are derived as 




2£'ik — J‘^{mcmf + E^) 


3J N f-' 2v^k - 


k,± 


Ik 


2mc 

~r 




2Eii^mf — J^mc{mcmf + E^), 


k,± 




2k 


1 

N 


K± 

k,± 

‘2mf 

J 


Ew^rric — ^mfirricmf + E^) — 2m 


c^ki 




2k 


( 10 ) 


where tanh(i?^/2J’). In order to draw the phase 

4ii/k 

diagram on the J-T plane, the free energies of the pure AFM 
phase and KP phase should also be calculated. 

The free energy of the KP phase is given by 


- ^^ln[2cosh(£;±/2T)], (11) 


with E^ = ^ 12 ± |ek|\/ek + and the 

equation for E: 

4 _ ^ ia.n\i[E^/2T) kkl 

^k + 

As T approaches Kondo temperature Tk, E —0, then E^ 
|efc|, E^ —> 0, therefore Tk of the KP phase is determined by 

4 _ 1 \" tanh(|ek|/2rR-) 

3J~2N^ kkl ■ ^ ^ 

For the pure AFM phase, the corresponding free energy den¬ 
sity is written as 

= Jmcmf - 2Tln[2cosh(£;+/2T)] 

2T 

-^V (13) 

k 

with the self-consistent equations for the AFM order parame¬ 
ters; 


At finite temperatures, the existence of thermal fluctuations 
may shift the parameter region of the CE phase. For simplicity 
and without loss of generality, we consider half-filling case 
ric = I with t' = 0. Since in this situation, the mean-field 
Hamiltonian has been diagonalized with the spectrums EE^, 


me 


Jrrif 
N 

mf — - tanh( Jmc/4T) = 0. 


tanh(i?j^ /2T) 

k ^ 
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When T approaches the Neel temperature T/v, mf,mc 0, 
then —>■ | Ck |, 0, thus the equation determining Tn 

is derived as 


'yiTjq 1 tanh(|ek|/2lAr) 

M ■ 

In order to derive the finite-temperature phase diagram, the 
free energies of the AFM phase, the KP phase and the CE 
phase are compared. In Fig. 1(c), the free energies of the three 
phases are plotted as functions of kondo coupling strength J 
at temperature T = 0.2t. It can be seen that the CE phase 
is stable in the coupling range Jci < J < Jc 2 , as it exhibits 
lowest free energy, and the AFM and KP phase occur in the 
coupling region J < Jd and J > Jc 2 , respectively. The 
phase boundaries Jd and Jc 2 now vary with temperature, and 
are crucial to determine the phase diagram on the J-T plane. 
Alternatively, we can consider the characteristic temperature 
separating CE phase with AFM phase as a function of J. On 
this boundary, V approaches zero, so it appears as the Kondo 
temperature in CE phase. Using Eq. (10), the equations 
determining are reduced to 

4 1 tanh(£;^/2r^) 2Eik - J’^rricmf _ ^ 

3J ^ 2^Et^ - E'i^ ^ ’ 

me rrif tanh{E^/2T^) _ 

T " Iv 4- 4^+ " ’ 

k 

2m/= tanh( Jmc/4T^), (15) 

where i^ik = + m'j), i? 2 k = 

Jmesje^ + E^ = + J^m^^ and 

i?k = '^"ic/2. 

On the boundary between CE phase and KP phase, the 
AFM order vanishes, so this phase boundary line corresponds 
to the Neel temperature T'pf. On this edge, m/ /me remains fi¬ 
nite. The self-consistent equations determining T'j^ with vary¬ 
ing J are simplified to 

tanh(i?^ kk| ^ 

4£;± + ’ 

2 1 ^ tanh(E;^/2rjV) ^ , 2Ei^j - 

J 4E;± 

E^^--iPV^l2-2tl ^ 
(16) 

where i5ik = Ck + !% E' 2 k = 12, E^ = 

71\/4 + kk|\/ek + and 7 = mflme- 

The critical lines Tk, Tat, and (which are all cal¬ 
culated as functions of J) necessary to determine the finite- 
temperature phase diagram at half-filling case are illustrated 
in Fig. 3. Both the Kondo temperature and Neel temperature 
show two distinct parts. In weak Kondo coupling region, the 


^ ^ taiil4E^/2^ 

J 4E;± ^ 



FIG. 3; (Color online) Finite temperature phase diagram at half¬ 
filling. In addition to the pure AFM phase and KP phase, a coex¬ 
isting phase emerges in the area between the lines and Tn 
and Tk are Neel temperature and Kondo temperature, respectively. 
Three phases converge to a single point on the J-T plane. 



FIG. 4: (Color online) Temperature dependence of staggered magne¬ 
tization M and KS strength V at fixed Kondo couplings. 


Neel temperature first increases from zero to a maximal value 
T/v = 0.238f at J = 2.85t, then diminishes rapidly down 
to zero at a critical Kondo exchange Jc 2 = 3.4161f. The re¬ 
duction of T/v in the intermediate coupling region is due to 
the spin fluctuations caused by the KS in the CE phase. For 
the Kondo temperature in the KP phase, it grows rapidly with 
increasing J from Tk = 0.238t at J = 2.85f , while inside 
the CE phase, it shows a steep reduction down to zero from 
J = 2.85f to Jei = 3.1498f. The CE phase exists in the nar¬ 
row region between and T^, which contracts continuously 
as the temperature rises, then finally disappears. 

Notably, the four lines intersect each other at J = 
2.85f,T = 0.238t, showing that the AFM phase, the KP 
phase and the CE phase converge to a same point on the 
J-T plane, similar to that reported in the case of far away 
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from half-filling, where the ferromagnetism and KS coex¬ 
ist [19, 20]. This feature of the phase diagram is similar to 
that derived by dynamic mean-field theory [2]. In fact, the CE 
phase deduced by previous mean-field treatments [6, 14] does 
not converge together with the AFM and KP phase to a single 
point in the phase diagram, in contrast, the CE phase dimin¬ 
ishes as temperature rises, and disappears inside the AEM and 
KP phase. Therefore, the existence of this convergent point 
deserves further verification beyond simple mean-field treat¬ 
ments. 

The staggered magnetization M and the Kondo hybridiza¬ 
tion parameter V are calculated as a function of J for given 
temperatures, illustrated in Fig. 2(b), and as functions of T 
but constant J, demonstrated in Fig. 4. On the edge between 
AFM phase and CE phase (denoted by Jd in Fig. 2 and by 
T'j^ in Fig. 4 (b)), M varies continuously showing a kink, then 
decreases and approaches zero on the upper edge Jc 2 or Tj^; 
while on the boundary between CE phase and KP phase (Jc 2 
in Fig. 2 and T'j^ in Fig. 4 (b)), the KS strength V also varies 
continuously with a kink and then decreases and disappears 
when approaching the lower edge {Jd or T^). The suppres¬ 
sion of AFM order and Kondo hybridization by each other in 
CE phase is owing to the competition between them. 


V. GROUND-STATE PHASE DIAGRAM CLOSE TO 
HALF-FILLING 


In above sections, the CE phase is studied in the half-filling 
case with only NN hoping t among conduction electrons, and 
the systems is in an insulating phase. While away from half¬ 
filling or with NNN hoping t' or beyond, the system no longer 
possesses particle-hole symmetry, thus in addition to the sin¬ 
glet hybridization 14 , the triplet hybridization V) between the 
conduction electrons and local moments plays an important 
role. Consequently, the system may possess enriched phase 
transitions and phase diagram. In contrast to the mean-field 
methods in earlier works [6, 14], the optimized mean-field de¬ 
coupling we employed here can be naturally generalized to in¬ 
clude the influence of Uc and t', and various phases and phase 
transitions between them can be discussed explicitly. 

In general case, the quasiparticle spectrums of CE phase 
have to be derived by diagonalizing "Hko- numerically, and the 
unitary transformation between the quasiparticles and c and / 
fermions can also obtained through this computation. The six 
self-consistent equations for CE phase are derived by fitting 
the number of c and / fermions to ric and 1, respectively, and 
by the definition of mean-filed parameters 14 , V], rric, ruf. 
Each equation is expressed in turns of the matrix elements of 
the unitary transformation. These equations are solved iter¬ 
atively until convergence is reached, then the energy of CE 
phase is obtained by summing the excitations below Fermi 
level. For pure AFM phase and KP phase, since the analytic 
spectrums exist, these two phase can be solved by minimiz¬ 
ing their ground-state energies. In the AFM phase, the con¬ 
duction electrons and /-fermions are decoupled, with A = 0, 
nif = 1/2, causing a smooth dispersions E^i = ±i Jmc of 
local spins. 


In order to determine the phase boundaries among the CE 
phase, AFM phase and KP phase, we have to develop an ef¬ 
ficient theory. Considering when J —Jc 2 , the parameters 
mf,mc, Vt ^ 0, the quasiparticle spectrums of CE phase can 
be expressed by two parts; one is the function of 14, A, /r, and 
the other can be perturbed in the first-order of to/, TOc, 14. To 
do this, we rewrite the mean-field Hamiltonian to the form 


= iVeo + ^'$L?fk.$k., (17) 

the operator are redefined as $kCT = 

(cko- ha Ck-fQo- /k+Qcr)^, and the Hamiltonian matrix 


ek-/r -jJVs ^Jrrifa jJVtcr \ 

^ -jJVs A iJVta -\JmcCr 

Urufa \^JVta Ck+q --jJVs 

\ jJVta -^JrricCr -|J14 A ) 

_ f A\^ (Ti?k \ 

crHk Alk+q j 

Using the Bogoliubov transformation 


( C]^a\ _rj ( Ckko- \ “k -fk \ ( Ckko- \ 

V ha J '^{ha J V «k Uk V /3ka J 

with Mk + ^k “ 1’ block diagonal parts Ak and Ak+q 
are diagonalized: C/j^AkUk = diag(£^^*^,= Ak, 

^iq^k+Q^k+Q ~ ’^^®'S('®k-rq-i-’-^k-rq-) = ^k+q^ 
where the dispersions are functions of 14, A, p, and are equal 
to the spectrums in the KP phase: 

4± = ^[ek-M + A±y/(ek - M - A)^ + QJ^UZ/d]. (18) 

To construct a global zero-order unitary transformation with 
Bogoliubov transformation 


Mk = 


C/k 0 \ 

0 C/k+q J ’ 


which acts on the Hamiltonian matrix leads to 




Ak crDk \ 
aD^ Ak+q J 


The off-diagonal elements hybrid the zero-order eigenstates 
with each other, hence bring a correction to the dispersions: 

Ek± = + El^l{mc,mf,Vt), where are 

easily obtained by these off-diagonal elements using perturba¬ 
tion theory. To the second order of (toc, to/, 14), the ground- 
state energy density near Jd can be divided into two parts 
E'^^ = Ef\Vs,\h} + Ejj^\mc,mf,Vt), where 

= I E + pVs - A + 

k,± 

= I E ^(-4±)4± - + Jrn^rnf. (19) 

' k,± 
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FIG. 5; (Color online) Energies of AFM phase KP phase 

and CE phase E'^^ vs Kondo exchange J. All energies are in unit 
of NN hopping t. 


EIG. 7: (Color online) Staggered magnetization M, Kondo hy¬ 
bridization Vs and triplet hybridization Vt as functions of Kondo 
coupling J. 


t'/t=0 t'/t=0.1 



n 

c 


n 
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EIG. 6: (Color online) Ground-state phase diagram of KLM near 
half-filling. As ric reduces, the collapse of KS at Jd and the mag¬ 
netic transition at Jc 2 get closer and finally converge. Thick and thin 
lines denote first and second order phase transitions, respectively. 



n 
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EIG. 8: (Color online) Staggered magnetization M, Kondo hy¬ 
bridization Vs and triplet hybridization Vt on Jd as functions of 
ric for CE phase. At ric = 1, Vs and Vt approach zero (see Fig. 
7a), while on the triple point (n*=0.8228 for t'ft = 0.1), M and Vt 
disappear. 


Minimization of Eg^'^ with respect to 14, A, p gives rise 
to three self-consistent equations, while differentiating Eg^'^ 
with {mc,rnf,Vt) produces other three equations. At Jc 2 , 
{nicmf^Vt) —>■ 0, but mf/rric and Vt/rric remain finite. 
Solving the six equations, Jc 2 and the value of Vs, X, p, 
m f Irric, Vt /rric on this boundary are calculated. 

In Fig. 5, the energies of the pure AFM phase, the CE phase 
and the KP phase are plotted with varying J, and the derived 
critical Kondo coupling Jd between AFM and CE phase on 
which the two phases have equal energy has been given as a 
function of ric and t' in Eig. 6. At half-filling, the pure AEM 
phase is separated with the CE phase by a second-order phase 
transition at Jd, and on this boundary, M varies continuously, 
while Vs and Vt approach zero, see Eigs. 7(a). Eor ric < 1, 


the transition at Jd changes to a first-order one, as indicated 
by the kink in Eg (Eig. 5(b)) and the discontinuity of M, 14 
and Vt at Jd (Eig. 7(b)). Note that Vs and Vt remain finite at 
Jd for ric < 1. At Jc 2 , a second-order phase transition be¬ 
tween CE and KP phase takes place, as seen by the tangency 
of ground state energy at Jc 2 in Pig. 5. Eg, Vs, Vt, M dW vary 
continuously with J at Jc 2 , at which Vt, M approach zero. 
Moreover, the NNN hoping t' can shift both boundaries. We 
find a sudden jump of Jd at ric = 1 (see Pig. 6), this fea¬ 
ture can be understood by the discontinuity of M and E^^ 
vs ric on Jd ( see Pig. 8) and may attribute to the change 
of topology, i.e., from the first order transition at ric < 1 to 
second-order transition at ric = 1- Only Jc 2 represents a real 
phase transition, because at this boundary the staggered mag¬ 
netization M vanishes from AEM to KP phase. 
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(a) J=2.5t 


(b) J=3t 


(c) J=3.5t 


tal observations for Celna and CeRhi_a;Coa;In 5 [12, 13], and 
YbRh 2 Si 2 under Co and Ir doping and external pressure [7,8]. 

The triple point (n*, J*) in our phase diagram can be 
shifted by t' jt, so in our mechanism this offset-to-onset transi¬ 
tion can also be generated by varying t'/t, similar to that pro¬ 
posed in Ref. [10]. Chemical or external pressure may simul¬ 
taneously change t'/t and ric, so which path cut in our phase 
diagram corresponding to these experiments is not clear. Ex¬ 
perimental studies of the existence of n* may be particularly 
interesting. The KP phase possesses larger Fermi surface than 
AFM phases, consequently, the transition from AFM to KP at 
Jc when ric < n* may induce a abrupt change of Hall coef¬ 
ficient as in YbRh 2 Si 2 , where the FSR was observed via Hall 
effect at the onset of magnetic transition [8]. 


VI. CONCLUSION 


FIG. 9: (Color online) Spectrums and Fermi surfaces of (a) AFM 
phase, (b) CE phase and (c) KP phase at ric = 0.95 and t'/t = 0.1. 
The shaded areas are the occupied Fermi sea, the white zones denote 
Fermi holes. 

The spectrums and Fermi surface topological structures of 
these three phases are shown in Fig. 9. The itinerant elec¬ 
trons and local moments are completely decoupled in pure 
AFM phase, leading to a hole-like Fermi surface around 
(7r/2,7r/2,) at t' = 0.1, which is consisted of the conduc¬ 
tion electrons only and the Futtinger volume equals UcSp 
, where Sp is the area of Brillouin zone; while in the CE 
phase, the hybridization of conduction electrons and local mo¬ 
ments constructs a hole-like Fermi surface around (0,0) and 
(tt, tt) points, with Fermi surface volume ncSp- The change 
of Fermi surface topology between AFM and CE phase indi¬ 
cates a first-order Fifshitz transition at Jci- For the KP phase, 
a hole-like Fermi surface exists around (tt, tt) with Fermi sur¬ 
face volume containing both c- and /-fermions. 

We have calculated Jd and Jc 2 as ric varies from 0.7 to 1. 
Notably, as ric decreases, Jd and Jc 2 approach each other and 
finally intersect at n* = 0.8228, J* = 2.2181f for t'/t = 0.1, 
and at n* = 0.8301, J* = 2.4604f for t' = 0, respectively. 
This result indicates that the CE phase region reduces with 
decreasing ric and finally disappears at a triple point. As ric 
decreases further, the AFM and KP phases are separated by a 
first-order transition at Jc, at which the AFM and KP phase 
shares equal energy. The ground state phase diagram of the 
KFM are summarized in Fig. 6 on the ric-J plane for given 
hopping parameter t' jt = 0 and 0.1. When ric > u,*, the 
collapse of 14 at Jd and the magnetic transition occurring at 
Jc 2 separate. While ric < u,*, the KS collapses precisely at 
the magnetic transition point. Similar phase diagrams are also 
given by mean-field treatments, Gutzwiller approximation and 
variational Monte Carlo approach [10, 15, 17]. This offset- 
to-onset transition between Kondo breakdown and magnetic 
transition as ric decreases may account for the experimen- 


In summary, we have performed an optimized mean-field 
decoupling of the Kondo lattice model near half-filling at both 
zero and finite temperatures. In addition to the pure AFM 
phase in weak Kondo coupling range and the Kondo paramag¬ 
netic phase in relatively strong coupling range, a distinct phase 
coexisting AFM order with Kondo hybridization arises in the 
intermediate Kondo exchange region, and the ground state 
phase diagram has been determined as function of the Kondo 
coupling, electron concentration and electron hopping. In par¬ 
ticular, for the coexisting phase, we found a finite staggered 
triplet hybridization between local moments and conduction 
electrons. We also develop an efficient method to calculate the 
phase boundaries. The characteristic parameters and Fermi 
surface structures of these phases and the phase transitions 
between them have been discussed explicitly. We have further 
found a mechanism explaining the offset-to-onset transition 
between Kondo breakdown and magnetic transition, which is 
driven by the decreasing of electron number ric. This mech¬ 
anism may account for the separation of the two transitions 
in YbRh 2 Si 2 under doping, and the existence of triple point 
(n*, J*) in the phase diagram deserves deep experimental in¬ 
vestigation. At half-filling limit, a finite-temperature phase 
diagram has been determined on J-T plane. As temperature 
rises, the region of this coexisting phase diminishes contin¬ 
uously then finally converges to a single point, together with 
the pure AFM phase and KP phase, which may require further 
theoretical and experimental verification. 
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